Scaling law of the plasma turbulence with non conservative fluxes 
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It is shown that in the presence of anisotropic kinetic dissipation existence of scale invariant power 
law spectrum of plasma turbulence is possible. Obtained scale invariant spectrum is not associated 
with the constant flux of any physical quantity. Application of the model to the high frequency part 
of the solar wind turbulence is discussed. 
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Dissipation range of incompressible hydro dynamic tur- 
bulence was extensively studied by different authors 
QEHSEH0, S H due to the fact that smallest 
scale perturbations display strong intermittency even at 
Reynolds numbers so low that there is no basis for fractal 
cascade. The kinetic energy spectrum E(k) of the hydro- 
dynamic turbulence in the far dissipation range behaves 
as 



E(k) ~ k ai exp [-a 2 (k/k d ) r 



(1) 



where kd is Kolmogorov dissipation wave number, oi and 
0:2 are some constants and 1 < n < 2 (see, e.g., Q and 
references therein). 

In the case of plasma turbulence existence of various ki- 
netic mechanisms of dissipation makes the situation much 
more complicated. For instance, observations of the so- 
lar wind turbulence 0, 0, strongly suggest steep 
power law spectrum of the magnetic field fluctuations for 
the frequencies higher then ion cyclotron frequency for 
which kinetic mechanisms of dissipation are dominant. 
In contrast with viscosity kinetic dissipation of plasma 
waves in the presence of background magnetic field is 
usually strongly anisotropic. In the presented paper we 
show that in the presence of: (i) anisotropic kinetic dis- 
sipation; and (ii) if the nonlinear transfer is governed 
by the scattering of the plasma waves by low frequency 
ones, then one should expect scale invariant power law 
spectrum of the plasma turbulence. It should be empha- 
sized that the scale invariant spectrum is not associated 
with the constant flux of any physical quantity due to 
the presence of kinetic dissipation. 

The general equation that governs the evolution of any 
averaged characteristic Z of the homogenous turbulence 
which is conserved by nonlinear interactions (in the case 
of hydrodynamic turbulence Z is usually associated with 
energy density) in the wave number space has the form 



d 4 = J + D + S, 



(2) 



where S and D describes the source and dissipation of 
Z and J accounts for the accumulation of Z due to the 
nonlinear interactions among the various wave number 
components of the turbulent field, such as velocity and 
magnetic fields. If nonlinear transfer term J serves only 
to redistribute Z and not changes total amount, then one 
can define the flux field F in the wave number space fl3| 



J=-V-F(k), 



(3) 



so the property of conservation is automatically fulfilled. 

Further we assume that nonlinear interactions are lo- 
cal in the sense that the most contribution to F(k) is 
from nearby regions of the k space. We make diffusion 
approximation to the wave number space transport, i.e., 
we assume that the flux can be presented as 



fi(k) 



-Di 



dQ_ 

' dkj ' 



(4) 



where are diffusion coefficients and Q is potential, 



that in general case is some function of k and Z. 
the inertial range Eq. J2J takes the form 



So in 



8Z_ 
dt 



_d_ 

dk t 



Dij(k,Z) 



dQ(k,Z) 



dk-j 



(5) 



The diffusion approximation for isotropic hydrod 
amic turbulence was first introduced by Leith 
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Afterwards the same concept was successfully applied 
to the plasma turbulence both in the strong [1J| and 
weak 1 151 turbulence regimes. For isotropic hydrody- 
namic turbulence Z corresponds to energy density £(k), 
Dij ~ k 9 / 2 5ij and Q ~ £ 3 ^ 2 13], where 5ij is kronecker 
delta. Combining Eqs. ©-(0} one can readily obtain the 
famous Kolmogorov spectrum £ ~ fc -11 / 3 for the inertial 
range of the hydrodynamic turbulence. In the case of 
the plasma turbulence the situation is more complicated. 
Existence of the background magnetic field usually leads 
to the anisotropy of the nonlinear cascade. 

Consider plasma turbulence in some frequency range 
where there is no source of the turbulence (S — 0) and 
there exist kinetic dissipation of plasma waves. As it was 
shown in Refs. 0, 01 , energy transfer in the dissipation 
range of hydrodynamic turbulence is dominated by non- 
local triads in which one leg is in the energy-containing 
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(low-k) range. Similarly, we suppose that the strongest 
nonlinear interaction is the scattering of the high fre- 
quency waves by low frequency ones from the inertial 
range of the plasma turbulence. In the weak turbulence 
theory this process requires fulfillment of the reso- 
nant conditions 



ki 



K, U)\ — Ll>2 + 0, 



(6) 



where kj. 2) ^1,2 and K, Q are wave numbers and fre- 
quencies of high and low frequency waves respectively. 
Due to the fact that uj ^> fl and fci,2 3> K, the change of 
the wave number is small in the unit act of the scattering 
( |ki — k2 1 <C |fei,2|) and therefore diffusion approximation 
is applicable (see, e.g., 0,H3). 

We incorporate the dissipation in the model as follows: 
we assume that the dissipation of the waves is negligible 
compared to the nonlinear interaction if the angle of the 
propagation with respect to the external magnetic field 8 
is less then some angle 80, whereas the opposite limiting 
case takes place for 8 > 9q. This model seems reason- 
able for the transverse waves in collisionless plasma, such 
as whistler waves and electromagnetic (ordinary and ex- 
traordinary) waves in electron-positron plasma 0, . 
Indeed, in this case the main mechanisms of dissipation 
are Landau and cyclotron damping. Both Landau and 
multiple cyclotron resonances do not affect the transverse 
waves for parallel propagation with respect to the back- 
ground magnetic field (see, e.g., 0]), whereas for rela- 
tively large angles of propagation both mechanisms can 
be on work. In the case of Landau damping this is caused 
by the fact that ambient propagating wave has nonzero 
electric field component parallel to the background mag- 
netic field - necessary condition of Landau damping. 

The second assumption of the presented model, that 
the nonlinear transfer is governed by the scattering of 
plasma waves by low frequency ones is also satisfied for 
both kinds of mentioned above wave modes. For of 
whistler waves in the solar wind, as it is shown below the 
strongest nonlinear process is scattering of whistler waves 
by low frequency magnetohydrodynamic waves from in- 
ertial range of the solar wind turbulence. Similarly, the 
strongest nonlinear process that governs evolution of elec- 
tromagnetic waves in electron positron plasma is their 
scattering by low frequency Langmuir waves (see e.g., 
[20| and references therein). 

It should be noted that presented model is not valid 
for Alfvenic turbulence. Alfven waves have low frequency 
compared to the ion cyclotron frequency and therefore 
they are not affected by cyclotron damping. Landau 
damping of Alfven waves also have unusual properties 
|2*H incompatible with the presented model. Addition- 
ally, the second assumption about nonlinear transfer is 
also violated in this case (see, e.g., [i^L 

Note, that if 8q is not extremely small, presented model 
implies that the wave with 8 — should take part in 
many scattering events before it can be transferred to the 



dissipation area. This circumstance allows us to use Eq. 
|JS| for the conical area in the wave number space 8 < 8 , 
and take the dissipation into account by requesting Q(k) 
to vanish at 8 = 8 . 

In the case of Alfvenic turbulence numerical simula- 
tions 23] as well as analysis of three and four wave 
resonant conditions provides that turbulent cascade is 
strongly anisotropic. In the weak turbulent regime there 
is no cascade in the direction parallel to the background 
magnetic field at all j^. When the dispersion of both 
high and low frequency waves can be considered as nearly 
isotropic for 8 < 8q, then one can expect different dif- 
fusion coefficients for the directions parallel and per- 
pendicular to the wave vector k of the high frequency 
waves. We consider scale invariant diffusion, i.e., assume 
Dij(k,Z) = dijSijk 01 *, and Q ~ k^Z^ = k 02 E(k)^ 
where dij are not functions of k and Z. 

With these assumptions Eq. (J5J reduces to the follow- 
ing one 



k 2 dk \ 11 dk 



D± d / . dQ 
k 2 sm8d8 V m dd 



= 0, (7) 



with boundary condition Q(k, 80) = 0. Here D\\ = d\\k a u 
and D± = dj_k a± k" 11 a± , and k s denotes the minimal 
wave number for which the formulated model is valid. Let 
us first consider the case when parallel and perpendicu- 
lar diffusion coefficients have the same scaling law, i.e., 
an = a± = a. Using standard methods of variable sepa- 
ration, i.e., representing potential as Q(k, 8) = R(k)^f(6) 
we obtain for the solution of Eq. J7J 



Q=J2 B m k- Cm P Vm (cose), 



(8) 



where coefficients B m are determined by the source of the 
turbulence at small wave numbers k = k s , P„ m (cos 8) are 
Legendre functions of the first kind and v m are the solu- 
tions of the eigenvalue problem P u (cos8q) = 0, arranged 
in order of increasing magnitude and 



-a - I 



'(a + 1) 2 + 4^ 



(9) 



If 6*o < tt/2 the value of the first eigenvalue v\ can be 
approximated as [2^ 



2.405 I 



(10) 



80 2 

For high wave numbers (k 3> k s ) the leading term is 

Q w Bxk- Cl P Ul (cos9). (11) 

Note that without kinetic dissipation isotropic solutions 
of Eq. is Q - k-"- 1 [E(k) - fc-(«+ 1 +^)/T2] 5 w hich 
correspond to constant flux of Z. In contrary, obtained 
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scale invariant asymptotic solution [see, Eq. iJHJ] is 
not associated with constant flux of any physical quan- 
tity, due to the presence of kinetic dissipation. Alterna- 
tively, in contrast with hydrodynamic turbulence, where 
the spectrum in the dissipation range is exponential, ob- 
tained result shows that the energy spectrum should de- 
crease as a power law in the dissipation range if the dif- 
fusion coefficients have the same scaling law. 

Now consider the case a± — a\\ = Aa / 0. Using the 
same technique of variable separation and introducing 
new variables K = k/k s and P(K) = k 1+a \\/ 2 R(K), Eq. 
yields 



K" 



: d 2 P 
&K 2 



an (aii + 2) 



v.. 



2 ^L^Ao 



P = 0. 



(12) 



If Aa < 0, for K 3> 1 one can drop the second term 
in the squire brackets. This yields the result that coin- 
cides with the result of the isotropic case Q ~ fc -Q " -1 , 
i.e., asymptotically dissipation has no influence on the 
cascade. 

On the other hand, when Aa > 0, for K 3> 1 one can 
drop the first term in the squire brackets. Obtained equa- 
tion can be solved in the terms of modified Bessel func- 
tions. Using asymptotic properties of modified Bessel 
functions 25] we obtain 



-(a||+l)/2-Aa/4 



exp 




-L^Aan/2 



(13) 

Consequently, in the case under consideration the spec- 
trum is exponential. Although it should be noted, that 
the decay is more soft than in the hydrodynamic turbu- 
lence when Aa is relatively small (Aa < 2). 

One of the possible application of the presented model 
is high frequency part of the solar wind turbulence spec- 
trum. Various spacecraft observations show the presence 
of persistent magnetic fluctuations in the solar wind over 
a broad range of frequencies 0,0, El- For low frequen- 
cies (/ < 10 -2 — 1CU 3 Hz) the magnetic field spectrum 
vary as approximately Bjvf(/) ~ / • F° r higher fre- 
quencies, up to proton cyclotron frequency (/ ~ 0.1 — 1 
Hz) the Kolmogorov spectrum /~ 5 / 3 is observed. This is 
believed to be the inertial interval of the solar wind tur- 
bulence. The change of slope and rapid decrease in the 
intensity near the ion cyclotron frequency is usually con- 
sidered to be due to the absorption of Alfven waves by ion 
cyclotron damping or Landau damping [2fij . At the fre- 
quencies, higher then ion cyclotron frequency, weak but 
persistent level of magnetic fluctuations, that can be well 
approximated by the power law spectrum f~ 3 is observed 
up to the electron cyclotron frequency. These fluctua- 
tions are usually associated with the whistler waves • 
The nature of this high frequency part of the spectrum 
remains unexplained. 

Whistler turbulence have been intensively studied by 
different authors both in strong (2^, l28l 123 ] and weak 



|30L l3l| turbulent regimes. If one assumes the exis- 
tence of the inertial interval of the whistler turbulence, 
then Kolmogorov- type dimensional analysis yields for the 
magnetic spectrum 28] EmQz) ~ fc~ 7 / 3 , that is incom- 
patible with observations (note that due to the relation 
EM(f)df ~ EM(k)dk, and taking into account Dopplcr 
shift and dispersion of whistler waves / ~ k 2 , observed 
/~ 3 spectrum corresponds to k~ v with v ~ 5 — 6 in the 
wave number space). 

There exist several different directions of the research 
for explanation of high frequency solar wind spectrum. 
The first approach [23, |32] is based on the fact that gov- 
erning equations of Hall magnetohydrodynamics besides 
energy, conserves two other second order (with respect 
to the field variables) quantities - magnetic and general- 
ized helicity |33j . Therefore, stationary Kolmogorov- type 
spectrum can be "driven" not only by energy cascade, but 
also by the cascade of magnetic and generalized helicitics 
|27|. In Ref. [34| short wavelength dispersive properties 
of the magnetosonic-whistler waves have been studied as 
a possibility reason of the spectrum steepening. Alterna- 
tive approach to the explanation of the high frequency 
magnetic fluctuations spectrum in the solar wind im- 
plies incorporation of the linear kinetic effects, such as 
Landau and cyclotron damping. It has been shown [sfij 
that simple incorporation of dissipation term to the en- 
ergy budget equation leads to sharp cut off of the energy 
spectrum. On the other hand, total ignore of dissipation 
leads to much more smooth spectrum compared to the 
observed one. 

The model considered in the presented paper could 
have important consequences for the explanation of the 
high frequency part of the solar wind spectrum. Whistler 
waves propagating along the background magnetic field 
are affected by neither Landau nor cyclotron damping 
on multiple ion cyclotron resonances. Based on the nu- 
merical solution of linear Vlasov equation the angle 
$o at which kinetic dissipation becomes dominant can be 
estimated as 9q ~ ir/6. The level of whistler wave fluc- 
tuations is low in the sense that (b^/B 2 <C 1, where 
(b^) is rms of the whistler wave magnetic field fluctua- 
tions and Bo is the background magnetic field. There- 
fore the study can be held in the framework of the weak 
turbulence theory. Possible nonlinear processes includes: 
(a) four wave resonant interactions of whistler waves [it 
can be shown that if (9o < 7r/3 then three wave reso- 
nances of whistler waves are absent, i.e., with this re- 
striction for all three waves resonant conditions similar 
to do not have nontrivial solutions] ; (b) induced scat- 
tering of whistler waves by ions; and (c) scattering of 
whistler waves by low frequency magnetohydrodynamic 
waves from the inertial range of the turbulence, i.e., three 
wave interactions which involves two whistlers and one 
magnetohydrodynamic wave. Detailed analysis of the 
nonlinear processes of the solar wind whistler waves will 
be presented elsewhere. Here we note that characteris- 
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tic time scales of these processes are respectively pro- 
portional to t q ~ A^ 2 , 77, ~ Af _1 , and r c ~ 1, where 
N(k) = £(k)/w(k) is the number density of whistler 
waves. Consequently, the strongest nonlinear process 
that should be responsible for the formation of the high 
frequency spectrum is the scattering of whistler waves by 
low frequency magnetohydrodynamic ones. This process 
conserves the total number of whistler waves jl5j |. and 
therefore Z = N(k). 

It can be shown that Alfven waves do not interact with 
whistlers through three wave resonances, whereas kinetic 
Alfven waves do [3(|. Another possibility is scattering 
of whistler waves by fast magnetosonic waves. Here we 
consider only the second possibility. 

Analytical calculations of diffusion coefficients are very 
complicated even in the incompressible limit. In the pre- 
sented paper we perform qualitative analysis of three 
wave interaction of whistler and fast magnetosonic waves 
and determine relations between diffusion coefficients 
that correspond to the observed spectrum. For this pur- 
poses we use equations of incompressible Hall magneto- 
hydrodynamics |33| 



V x B) x B] 



(14) 



d(B + V x V) 
dt 



-V x [(B - V x V) x V] , (15) 



where time and space variables are measured in units 
of ion giroperiod uj^} and ion skin depth Aj, respec- 
tively. Analysis of these equations shows that for triad 
interaction of whistler and fast magnetosonic waves the 
strongest nonlinear term is the second one on the right 
side of Eq. i|14(l . Taking this into account and noting 
that the wave number of whistler waves is much greater 
then the wave number of fast magnetosonic waves, then 
Fourier transform of Eq. (|14fl yields 



dt 



-fe||-B (k x K)=ik x [b / x (V x b™)] kJ (16) 



where index k denoted the Fourier transform and su- 
perscripts w and / indicates that corresponding values 
correspond to whistler and fast magnetosonic waves, re- 
spectively. 

Using helicity decomposition [31j one can apply the 
standard technique of weak turbulence theory developed 
for two types of interacting waves 0. But we do not 
perform this analysis here due to the reason that only 
thing that we need for further analysis is scaling index 
of the matrix element of interaction T that immediately 
follows from Eq. I|16l) : T ~ k 2 . In the framework of weak 
turbulence theory the dynamics is totally determined by 
linear dispersions of waves and matrix element of interac- 
tion. The method of finding the scaling index of diffusion 
coefficient was developed in Ref. which in the case 



under consideration yields a = 
Therefore, Eq. ijTT|l yields 

AT(k) - k- 



-1 and /3i = 0, 71 = 1. 



(17) 



where c\ — i/i(d^/d||) 1//2 . 

To obtain spectral index 5 of corresponding energy 
spectrum E(f) ~ f~ s , we note that 



E(f)~E(k)^~£(k)k^~N(k)k3. 



(18) 



Taking also into account that for whistler waves / ~ fc 2 , 
Eqs. <H7J) and (JTSJl yield 



ci-3 



(19) 



As it was mentioned above, according to observations 
S w 3. Taking also O = tt/6, Eqs. JTJ3 and JTSJl yield 
d±/d\\ ~ 5. Obtained result seems reasonable, due to 
the fact that in the magnetized media perpendicular cas- 
cade rate usually significantly exceeds parallel cascade 
rate (see, e.g., |38| and references therein). 

In the presented paper plasma turbulence in the pres- 
ence of anisotropic kinetic dissipation is considered. It is 
shown that if the nonlinear transfer is governed by the 
scattering of the plasma waves by low frequency ones, 
then development of asymptotic scale invariant power law 
spectrum of the plasma turbulence is possible. Obtained 
scale invariant spectrum is not associated with the con- 
stant flux of any physical quantity due to the presence 
of kinetic dissipation. Corresponding spectral index is 
given by Eq. © with m = 1 . Possible application of the 
presented model to the high frequency part of the solar 
wind spectrum has been analyzed. 

Author is immensely grateful to Alexander Bershadskii 
and George Machabcli for helpful discussions and sugges- 
tions. 
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